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ABSTRACT 

Doppler  radar  tomography  is  a  technique  which  uses  only  the  Doppler  profiles 
of  a  rotating  target  as  projections  in  two-dimensional  radar  tomography  to 
obtain  an  image  of  the  target.  The  work  demonstrates,  for  the  first  time,  the 
application  of  the  fractional  Fourier  transform  and  the  S-Method  in  improving 
the  image  resolution.  The  performance  is  shown  to  be  significantly  better  than 
that  available  with  traditional  Fourier  transform-based  methods. 
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Application  of  the  Fractional  Fourier  Transform  and 
S-Method  in  Doppler  Radar  Tomography 

Executive  Summary 

Radar  tomography  is  the  process  of  constructing  a  two-dimensional  image  from  one¬ 
dimensional  projections  obtained  from  radar  returns.  The  backscatter  from  radar  pro¬ 
vides  Doppler  information  of  the  scatterer  giving  the  frequency  characteristics  over  time, 
or  in  the  case  of  rotating  objects,  over  varying  aspect  angle.  These  frequency  domain 
projections  are  used  for  Doppler  radar  tomographic  imaging. 

Doppler  radar  tomography,  in  conjunction  with  other  signal  processing  techniques, 
can  also  be  used  to  estimate  the  physical  parameters  of  a  rotating  object.  For  helicopter 
targets,  the  characteristics  of  the  rotor  blades  (blade  length,  number  of  blades  and  rotation 
rate)  and  rotor  hub  features,  can  be  used  as  discriminates  for  automatic  target  recognition. 
The  ability  to  accurately  estimate  these  features  effects  the  ability  to  correctly  classify 
targets  quickly  and  is  directly  related  to  the  frequency  resolution  of  the  system. 

The  resolution  of  a  tomographic  image  depends  on  the  resolution  of  the  input  pro¬ 
jections.  Here,  we  strive  for  high  frequency  resolution  so  that  closely  spaced  scatterers 
can  be  resolved  and  hence  appear  in  the  tomographic  image.  However,  the  long  coherent 
processing  interval  needed  to  increase  frequency  resolution  introduces  a  frequency  spread 
(chirp)  for  each  scatterer  (due  to  the  scatterer  motion),  giving  ‘smeared’  projections  and 
hence  smeared  images. 

The  fractional  Fourier  transform  is  a  technique  that  can  be  used  to  identify  chirps  in 
a  signal  so  that  they  can  be  spectrally  compressed  using  the  S-Method.  The  result  is  high 
resolution  projections,  without  the  frequency  spread,  and  a  clear  and  sharp  tomographic 
image.  Long  integration  times  allow  closely  spaced  scatterers  to  be  easily  resolvable  and 
performance  is  improved  in  the  presence  of  noise  when  compared  to  the  traditional  Fourier 
transform-based  method,  which  is  confined  to  short  coherent  processing  intervals. 

This  report  successfully  demonstrates  the  application  of  the  fractional  Fourier  trans¬ 
form  and  the  S-Method  techniques  on  simulated  rotating  point  scatterers  and  real  heli¬ 
copter  data.  The  improved  resolution  of  the  projections  and  clarity  of  tomographic  images 
makes  identifying  the  position  of  scatterers  easier  and,  for  the  case  of  real  helicopter  data, 
the  blade  tip  positions  more  discernible. 

The  accurate  estimation  of  target  parameters  using  high  resolution  imaging,  is  expected 
to  improve  automatic  target  recognition  performance  and  can  be  used  on  a  variety  of 
targets  such  as  helicopter  rotor  blades,  propellers  on  propeller  driven  aircraft  or  slower 
rotating  objects  such  as  rotating  antenna  structures  on  ships. 
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1  Introduction 

Tomography  is  a  mathematical  technique  of  reconstructing  a  2-dinrensional  function  from 
1-dinrensional  projections  of  the  function  at  different  aspect  angles.  This  involves  the 
Radon  transform,  first  described  by  Johann  Radon  [1],  and  back-projection  algorithms 
to  reconstruct  the  original  function.  It  can  also  be  described  as  a  method  of  computing 
the  inverse  2-dinrensional  Fourier  transform  using  a  collection  of  1-dinrensional  projections 
which  can  be  measured  directly  from  an  object  of  interest.  In  X-ray  tomography,  which 
is  usually  referred  to  as  Computer  Aided  Tomography  (CAT),  the  projections  are  density 
measurements.  In  radar  tomography,  they  are  range  profiles. 

For  radars  of  high  bandwidth,  the  range  profiles  can  be  down-range  profiles  which  are 
basically  1-dinrensional  images  of  a  target  in  the  down-range  dimension.  For  narrow-band 
radars,  the  range  profiles  are  necessarily  cross-range  profiles  which  can  be  obtained  through 
coherent  Doppler  processing,  which  means  the  dwell  time  requirement  for  narrow-band 
must  be  greater  than  in  the  case  of  wide-band  radars.  The  trade-off  between  bandwidth 
and  time  is  a  manifestation  of  the  well-known  and  fundamental  time-bandwidth  product 
which  has  a  significant  role  in  setting  the  performance  limits  of  any  radar  system. 

The  theory  of  tomography  itself  has  been  well  described  by  many  authors,  such  as  in 
[2] .  Our  present  work  is  focused  on  the  use  of  Doppler  tomography  in  narrow-band  radars 
and  applying  time-frequency  processing  techniques  to  achieve  high-resolution  images  of 
rotating  targets.  D.  Mensa  [3]  was  among  the  first  authors  to  have  used  the  term  ‘Doppler 
tomography’.  Mensa  described  the  basic  ideas  of  Doppler  tomography:  all  scattering 
centres  at  the  same  cross  range  of  an  illuminated  rigid  target  impart  the  same  Doppler 
frequency  components  in  the  returned  signal,  and  hence  a  Doppler  profile  can  be  used 
to  represent  a  cross-range  profile.  As  the  target  rotates,  a  continuum  of  such  projections 
can  then  be  obtained  and  inverted  by  tomography  to  produce  an  image.  G.  Fliss  [4]  later 
applied  the  idea  to  the  imaging  problem  of  helicopter  rotor  blades  and  hubs,  using  a  simple 
single  frequency  continuous-wave  (CW)  instrumentation  radar.  Fliss  demonstrated  that 
it  is  possible,  using  Doppler  tomography,  to  image  the  tips  of  the  blades  which  are  the 
main  scattering  centers  of  the  blade  rotors. 

Resolution  of  a  Doppler  tomographic  image  directly  depends  on  the  resolution  of  the 
projections  used  in  the  tomographic  processing.  As  a  target  rotates,  the  instantaneous 
cross  range  position  of  a  rotating  scattering  point,  relative  to  the  center  of  rotation,  defines 
its  instantaneous  Doppler  frequency,  while  the  finite  coherent  processing  interval  (CPI) 
required  to  produce  a  projection  undesirably  introduces  a  finite  frequency  spread  for  each 
point  scatterer.  In  other  words,  each  point  scatterer  in  a  projection  is  represented  by 
a  chirped  signal  of  finite  spectral  extent,  rather  than  a  ‘pure  tone’,  causing  a  smearing 
in  a  projection,  which  in  turn  leads  to  a  degradation  of  the  tomographic  image,  if  not 
appropriately  processed.  The  longer  the  CPI,  the  finer  the  Doppler  filter  width  would 
become,  but  the  smearing  due  to  cross-range  motion  also  gets  worse. 

Known  techniques  for  spectral  compression  for  chirped  signals  include  bilinear  tech¬ 
niques  such  as  the  Wigner-Ville  Distribution  (WVD),  the  Cohen’s  class,  and  the  time- 
frequency  distribution  series  (TFDS)  as  discussed  in  V.  Chen’s  book  [5].  The  spectral 
compression  performance  of  the  WVD  is  excellent  except  for  the  existence  of  cross  terms. 
A  recently  proposed  variant  of  the  WVD  is  the  S-Method  [6] ,  which  could  remove  the  cross 
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terms  while  also  substantially  reducing  the  computational  cost.  Unfortunately,  when  the 
instantaneous  component  frequencies  are  closely  spaced,  cross  terms  are  still  unavoidable; 
and  this  is  indeed  the  case  in  Doppler  tomography.  Scattering  centers  may  line  up  with 
the  radar  line  of  sight,  causing  spectral  overlapping  and  crossing  of  the  chirp  components. 

Fortunately,  the  recently  developed  technique  of  fractional  Fourier  transform  (FrFT) 
has  become  available,  which  can  turn  the  spread-out  spectrum  of  a  chirp  signal  into  a 
well-compressed  spike  in  the  ‘fractional  frequency’  domain.  We  have  exploited  the  useful 
property  of  the  FrFT  and,  in  combination  with  the  use  of  the  S-Method,  developed  an 
algorithm  for  extracting  and  spectrally  compressing  multiple  chirp  signals,  to  produce 
clean  and  sharpened  projections  for  tomographic  processing.  The  application  of  the  FrFT 
to  Doppler  tomography  and  the  new  algorithm  of  extracting  and  compressing  multiple 
chirp  signals  are  the  main  contribution  of  this  work,  which  are  demonstrated  using  both 
simulated  and  real  data. 

An  outline  of  the  Report  is  as  follows.  Section  2  gives  an  overview  of  the  theory  of  2D 
tomography  in  general  and  the  system  requirements  for  Doppler  tomography  in  particular. 
Section  3  briefly  describes  the  fractional  Fourier  transform  and  its  applicability  to  Doppler 
tomography.  Section  4  describes  the  S-Method  and  why  it  is  useful  for  the  current  problem. 
Section  5  presents  the  algorithm  of  multiple  chirp  extraction  and  compression.  Section  6 
describes  and  discusses  the  results  when  applied  to  both  simulated  and  real  data.  The 
final  Section  contains  the  concluding  remarks. 


2  An  Overview  of  2D  Radar  Tomography 

This  section  summarises  the  theory  of  two-dimensional  tomographic  imaging,  explains 
what  is  meant  by  ‘Doppler  radar  tomography’,  and  finally  discusses  resolution  and  system 
requirement  issues  when  the  projections  used  in  the  tomographic  processing  are  based  on 
the  traditional  Fourier  transform. 


2.1  The  Basic  Theory 


The  essence  of  tomographic  imaging  is  the  Back-Projection  Algorithm  [2,  7,  8,  9,  10], 
which  can  be  formulated  as  a  general  2D  Fourier  analysis  technique  as  follows.  For  a  2D 
real-valued  function  q(x,  y),  as  shown  in  Figure  1,  the  projection  pg(r)  along  a  line  through 
the  origin  at  angle  6,  which  is  also  called  the  Radon  Transform,  can  be  defined  as 

/OO  coo 

/  q(x,y)  6(x  cos  6  +  y  sin  0  —  r)  dx  dy.  (1) 

-OO  J  —  OO 


The  spatial  Fourier  transform  Pg(fr )  of  pg(r),  with  spatial  frequency  in  the  r-  dimension 
denoted  as  fr,  is 


COO 

/  po(r)e~t2nf'rdr 

(2) 

J  r=— oo 
coo  coo 

/  /  q{x,y)  e~*2,r(/*  x+/i/ s/)  dxdy 

J — oo  J  —  OO 

(3) 

Qe(fxJy), 

(4) 
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A  projection 


Figure  1:  Projection  of  a  2D  function,  q(x,y). 


upon  the  substitution  of  (1).  Here, 


fx  =  fr  cos  8,  fy  =  fr  sin  8 


are  the  spatial  frequency  components  in  the  x-  and  y-  dimensions  along  the  slice  at  angle 
8,  and  Qg(fx,fy)  is  the  2D  Fourier  transform  of  q(x,y)  along  this  direction,  i.e.  a  slice 
of  Qoifx-  fy)-  The  spatial  frequency  fr  is  related  to  the  temporal  frequency  /  through 
fr  =  2 f/c.  The  Fourier  transform  of  a  projection  of  a  two-dimensional  function  is  thus  a 
slice  of  the  2D  Fourier  transform  of  that  function.  This  is  the  so-called  Projection-Slice 
Theorem. 


Once  Qo(fx,  fv)  is  constructed  from  many  such  slices  at  a  range  of  aspect  angle  8, 
q{x,  y )  can  be  obtained  from  the  inverse  2D  transform.  However,  noting  that  dfx  dfy  = 
fr  dfr  d8,  the  inverse  2D  transform  can  be  done  on  Qg(fr )  in  polar  coordinates,  as  follows: 


POO  POO 

q(x,y)  =  /  /  Qe{fx,fy)j2<xi*+yMdfxdfy 

J  —  OO  j  —OO 


'  — oo  J  — OO 
p  2n  poo 


/0  JO 


Qe(fr)  e^M-oose+ysinO)  ^  ^  ^ 
\fr\Pe(fr)  e^Mzcosfl+ysinf?)  ^  dQ 


’  0  LJ  —  oo 
—  1 


r  F-1  {\fr\  Pe(fr)}(xcosg+ysine)  d8, 

J  0 


T  1  (sgn(/r)  f.r  Pe{fr)}{xcosd+ys-me)  dd- 


(5) 

(6) 

(7) 

(8) 

(9) 
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The  subscript  (xcosd  +  ysin#)  means  the  inverse  Fourier  transform  (FT),  denoted  by 
(J--1),  is  calculated  at  r  =  xcosO  +  ysin.9.  This  J r~1  operation  is  where  most  of  the 
complexity  of  implementation  of  the  above  result  lies. 

There  are  two  known  approaches  to  calculate  the  above  integral  [2] .  The  first  approach, 
called  the  Filter  Back-Projection  method,  calculates  in  the  inverse  FT  in  (8)  directly 
(and  numerically).  The  second  approach,  called  the  Convolution  Back-Projection  method, 
applies  the  convolution  theorem  to  the  inverse  FT  in  (9).  In  general,  if  p(x)  is  any  real¬ 
valued  bounded  function,  then 

e-®/.  ix = i2*/  P(j). 

From  theory  on  Hilbert  transform, 

{sgn (/)}  =  h{x)  =  -  — , 
iirx 

then  the  convolution  theorem  leads  to  [2] 


J-  {sgn  (/r  )  fr  Pd{fr  )  }  (x  cos  9+y  sin  9) 


1 

i2ir 

1 

27 r2 


dpejr ) 
dr 


}  r=x  cos  0+y  sin  9 


dpejr ) 
dr 


r=x  cos  0-\-y  sin  6 


(10) 

(ii) 


i 

27 r2 

What  remains  is  to  integrate  this  distribution  (or  grid)  in  the  ( x ,  y )  plane  over  a  range 
of  angles  9  between  0  and  7r.  In  our  work,  we  have  used  the  Filter  Back-Projection  method. 


_ 1 _ dpojt) 

x  cos  9  +  y  sin  9  —  t  dt 


dt  (12) 


2.2  Doppler  Radar  Tomography 

In  applications,  the  projection  pg{r)  in  the  basic  theory  described  above,  can  be  obtained 
in  at  least  two  different  forms.  It  may  be  in  the  form  of  a  range  profile  of  a  rotating  target 
from  radar  with  a  high  sampling  rate,  which  may  be  the  case  for  radars  with  an  ISAR 
capability.  Alternatively,  it  may  be  in  the  form  of  a  Doppler  profile  of  a  rotating  target, 
which  can  be  obtained  by  a  narrow-band  Doppler  radar. 

The  technique  of  radar  tomography  which  uses  the  Doppler  profiles  of  a  rotating  target 
as  projections  may  be  called  ‘Doppler  radar  tomography’  and  is  the  main  subject  matter 
of  the  present  report. 

As  a  rigid  object  rotates  about  an  origin  such  as  that  depicted  in  Figure  1,  all  points  at 
constant  cross-range  (in  the  x  dimension)  offset  from  the  origin  will  have  the  same  Doppler 
shift,  assuming  the  receiving  antenna  is  located  in  the  negative  y  direction.  This  means 
that  the  magnitude  spectrum  of  the  Fourier  transform  of  a  sufficiently  long  segment  of 
samples  can  be  interpreted  as  a  cross-range  projection  of  the  object.  In  other  words,  if 
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s(t]tk );  tk  <  t  <  tk  +  Tcpi  denotes  the  received  time  signal  for  one  CPI  with  the  corre¬ 
sponding  averaged  angle  9}.  (say),  and  if  )  denotes  its  short-time  Fourier  spectrum 

at  slow  time  tk, 

ftk+Tcpi 

S{f,tk)=  /  w(t-tk)s(t,tk)e~z2nftdt,  (13) 

Jtk 

then  a  cross- range  projection  pgk  (r)  can  be  obtained  as 


Pek(r )  =  | S  (f,tk) | , 


(14) 


where 


2 LO  COS  Ok 
f  =  — ~r 


(15) 


and  to  is  the  rotation  rate,  still  to  be  estimated.  A  sufficient  number  of  such  projections 
at  different  angles  are  required  to  obtain  an  image. 


The  CPIs  are  not  required  to  be  non-overlapping.  In  fact,  numercial  experimentation 
has  shown  that  a  certain  degree  of  overlapping  between  adjacent  CPIs,  here  described  by 
a  factor  called  the  overlapping  factor  F,  can  lead  to  a  smoother  tomographic  image. 


The  estimation  for  the  rotation  rate  to,  which  is  usually  an  unknown  and  constant 
parameter  of  the  system,  can  be  achieved  by  a  number  of  different  methods  of  which  an 
iterative  procedure  based  on  some  quality  measure  of  the  resulting  image  is  one  possibility. 
Also,  object  translational  motion  compensation  must  be  done  prior  to  the  back-projection 
algorithm  processing,  which  can  be  implemented  simply  by  shifting  the  spectrum  S(f,tk) 
to  place  the  body  line  (which  is  also  usually  the  strongest  line)  at  zero  Doppler  before 
scaling  it  to  form  pgk  (r) . 


2.3  Image  Resolution  and  System  Requirements 

There  are  a  number  of  factors  determining  the  achievable  image  resolution  and  the  required 
system  parameters.  These  include  the  fundamental  limit  due  to  available  bandwidth 
expressed  by  the  point  spread  function,  and  various  factors  involved  in  the  formation  of 
projections.  These  factors  include  the  finite  sampling  rate  of  the  radar  system,  the  CPI 
used  for  a  projection,  as  well  as  the  rotation  speed  and  cross-range  dimension  of  the  object. 


2.3.1  The  Point  Spread  Function 


In  reality,  the  integration  limits  and  frequency  extent  are  all  necessarily  finite,  and  the 
available  ‘angular  bandwidth’  may  be  less  than  the  full  circle.  Equation  (7)  of  Section  2.1 
must  then  be  replaced  by 

f*7 r  I"  roc 

a(x,y)=  /  We(fr)\fr\Pe(fr)ei2^™se+ysine)dfr  ^  (16) 

Jo  U-oc 

where  Wq  ( fr )  defines  the  available  bandwidth 


We(fr) 


1,  ( fr,0 )  £  D  (some  domain) 

0,  otherwise 


(17) 
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and  sets  the  fundamental  limit  on  the  resolution  of  the  system.  The  2D  inverse  Fourier 
transform  of  Wg(fr)  is  called  the  point  spread  function  (PSF),  where  the  PSF  becomes 
more  focused  as  the  domain  D  gets  larger,  resulting  in  improved  image  resolution. 


2.3.2  Sampling  Rate  Requirement 

The  sampling  rate  (or  pulse  repetition  frequency  (PRF)  in  narrow-band  pulse  Doppler 
radars  which  collects  only  one  sample  per  transmitted  pulse)  needs  to  be  sufficient  to 
cover  the  entire  Doppler  frequency  extent  of  the  target.  More  exactly,  assuming  that  the 
body  Doppler  line  of  the  object  has  been  shifted  to  zero  frequency,  the  Doppler  frequency 
of  a  point  scatterer  at  the  largest  radius,  denoted  R-tip,  from  the  center  of  rotation  (i.e. 
the  ‘tip’  Doppler)  is 

f  _  ^vtrp  a  _  ZRt.ipU  Q 
Jt.ip  —  ^  COS  U  —  COSl/, 

which  takes  maximum  values  when  cos#  =  ±1  at  #  =  rmr;  m  =  0,1,2,....  By  the 
Nyquist  criterion,  the  minimum  sampling  rate  to  capture  the  highest  Doppler  frequency 
components  without  aliasing  (or  ambiguity)  is  simply 

PRFrmn  >  ■  (18) 

This  requirement  on  sampling  rate  applies  equally  to  all  types  of  Doppler  radar  tomogra¬ 
phy. 


2.3.3  CPI  Requirement  for  Fourier  Transform  Based  Systems 


In  Doppler  radar  tomography  based  on  the  Fourier  transform,  especially  the  Fast  Fourier 
Transform  (FFT),  the  length  of  a  coherent  processing  interval  (CPI)  is  limited  by  the 
condition  that  the  rotational  motion  of  a  point  scatterer  on  the  object  must  not  induce 
a  Doppler  spread  large  enough  to  spread  more  than  one  frequency  filter  width.  This  is 
because  the  frequency  filters  translate  directly  to  cross- range  bins,  and  any  spread  of  the 
signal  over  more  than  one  filter  would  translate  into  cross-range  ‘walk-off’,  or  smearing, 
in  the  resulting  image. 


Referring  to  Figure  1,  cross  range  is  given  by  x  =  r cos#,  hence  the  differential  change 
in  x  is 


dx  =  —  ru;  sin#  eft. 


Maximum  cross  range  walk-off  occurs  near  #  =  mr/2,  with  n  being  an  odd  integer.  If 
dt  represents  a  CPI  time,  now  denoted  by  Tcpi,  then  the  requirement  for  no  cross  range 
walk-off  during  a  CPI  translates  to  \dx\  =  ru>Tcpi  to  be  within  a  cross  range  bin,  for  all 
points  along  the  radial  direction  r. 

With  Rtip  again  denoting  a  largest  radial  distance  of  any  scatterer  on  the  object  relative 
to  the  centre  of  rotation,  Ns  =  Tqpi  x  PRF  being  the  number  of  samples  used,  which  is 
for  now  assumed  to  equal  to  the  FFT  length,  then  for  the  point  scatterer  center  at  Rup, 
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it  is  straightforward  to  show  that  the  requirement  becomes 


Ns  < 


PRF 


LO 


2  R, 


tip 


1/2 


=  PRF 


A 


\  1/2 

2  id  Vfip  ) 


{ 2PRF\  1/2  f  PRF \  1/2 


(19) 


as  a  condition  on  integration  time  Tqpi,  or  equivalently  the  number  of  samples  Ns  used  for 
integration,  for  smear-free  projections.  If  Ns  exceeds  this  limit,  resolution  in  a  projection 
may  be  improved  toward  the  tips  but  smearing  will  result  toward  its  center. 

Also  from  (19),  a  smear-free  CPI  time  for  each  projection  and  the  corresponding  angle 
of  rotation  of  the  scatterer  is 


Tcpi 
A  6 


Ns  PRI  = 
to  x  Tcpi  - 


1 

LO 


'  A  \1/2 
V  2  Rtip  ) 

\  y/2 

2  Rtip) 


(rad), 


(20) 

(21) 


which  is  independent  of  PRF,  but  dependent  on  radar  wavelength  A,  rotation  rate  to  and 
maximum  radial  extent  Rtip-  With  an  X-band  wavelength  and  a  radial  extent  between  4 
to  10  meters,  (21)  gives  Ad  between  2.2  to  3.5  degrees  for  each  projection. 

For  the  case  of  absolutely  minimum  PRF  given  by  (18),  the  result  in  (19)  gives  the 
required  maximum  number  of  samples  for  no  cross  range  smearing: 


Ns  < 


XRitp 

A 


1/2 

=  Npr , 


(22) 


which  is  an  upper  limit  independent  of  the  rotation  rate.  The  quantity  Npr  is  also  the 
maximum  number  of  points  spanning  the  cross  range  extent  of  the  target,  free  of  smearing, 
if  the  chosen  PRF  is  not  the  minimum.  It  also  means  that  a  fundamental  cross  range 
resolution  limit,  for  Fourier  transform  based  techniques,  is 

j  2  Rtip  f  ^Rtip\  ^  frtoX 

dx^  =  N^  =  \^)  -  (23) 

Note  that  this  resolution  limit  is  independent  of  to  and  PRF. 


2.4  An  Example 

For  a  simple  example,  two  rotating  point  scatterers  are  positioned  at  (—4,  —4)  and  (1, 1) 
in  the  (x,y)  plane  respectively.  They  rotate  at  a  rate  of  300  RPM  about  the  origin  and 
have  the  same  RCS.  The  following  radar  parameters  used  include  a  wavelength  A  =  0.03m, 
PRF  =  26  kHz,  and  a  dwell  time  of  approximately  0.2  sec  (one  cycle  of  rotation). 

Figure  2(a)  shows  the  smear-free  projections,  with  no  overlapping  between  adjacent 
CPIs  (i.e.  F  =  0),  of  two  point  scatterers  rotating  about  a  coordinate  origin,  with  the 
resulting  tomographic  image  shown  in  Figure  2(b).  While  more  discussion  will  be  given 
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(a)  Smear-free  projections  of  two  point  scatterers  with  F  =  0. 
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(b)  Tomographic  image  of  two  point  scatterers  using  the  smear-free 
projections 

Figure  2:  Smear-free  projections  and  resulting  tomographic  image,  based  on  the  Fourier 
transform.  Example  for  two  rotating  point  scatterers  rotating  at  300  RPM,  when  viewed 
by  an  X-band  radar  (X  =  0.03  m),  using  41  samples  over  a  0.2  s  period. 
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in  the  next  section,  it  suffices  here  to  indicate  that,  under  the  snrear-free  condition,  each 
CPI  consists  of  only  41  samples  and  the  pixel  size  of  this  image  is  approximately  0.29  m 
by  0.29  m. 

If  more  samples  are  used  in  a  CPI  in  an  attempt  to  increase  resolution,  smearing  in 
the  projections  and  the  resulting  image  occurs,  as  can  be  seen  in  Figure  3.  Here,  the  CPI 
length  is  increased  from  the  snrear-free  case  by  eight  times,  while  an  overlapping  factor 
of  0.6  has  been  used  to  obtain  a  reasonable  number  of  projections  for  imaging.  Figure 
3(a)  indicates  that  the  frequency  spread  is  most  pronounced  around  zero  Hertz  (when 
the  scatterer  is  moving  in  the  pure  cross-range  direction  relative  to  the  radar)  and  the 
least  at  the  tip  values  (when  the  scatterer’s  velocity  is  along  the  line-of-sight).  It  should 
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(a)  Smeared  projections  of  the  two  point  scatterers,  and  F  =  0.6. 


Cross  Range  (m) 


(b)  Tomographic  image  of  two  point  scatterers  using  the  smeared  pro¬ 
jections 


Figure  3:  Smeared  projections  and  resulting  tomographic  image.  The  same  scatterers  as 
used  in  Figure  2,  using  328  (41  x  8 )  samples  over  a  0.2  s  period. 


also  be  pointed  out  that  ‘smearing’  in  Figure  3(b)  means  the  ‘blobs’  representing  the  point 
scatterers  spread  over  multiple  pixels  in  that  figure,  compared  to  the  snrear-free  case,  whose 
overall  blob  size  is  smaller.  In  fact,  smearing  occurs  mostly  in  the  direction  tangential  to 
circular  path  of  rotation.  This  is  characteristic  of  the  tomographic  processing  and  partly 
due  to  the  fact  that  the  tip  scatterers  are  still  very  well  focused  in  the  projections. 


It  will  be  shown  in  the  next  sections  how  the  FrFT  and  the  S-Method  can  help  overcome 
this  limitation. 
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3  The  FrFT  and  its  Application 

This  relatively  new  tool  was  first  proposed  by  Narnias  [20]  in  1980  in  the  field  of  quantum 
mechanics,  which  has  been  discussed  extensively  by  many  authors  for  applications  in  signal 
processing  only  the  last  fifteen  years  or  so  [13,  14,  15,  16,  17,  18].  As  a  generalisation  of 
the  traditional  FT  in  a  rotated  time-frequency  plane,  it  has  great  potentials  in  the  analysis 
of  chirp  signals  which  become  ‘pure  tone’  signals  in  appropriately  rotated  time-frequency 
frames.  It  has  been  used  to  analyse  and  extract  multiple  linear  frequency  modulated 
(LFM)  audio  signals  from  bats  by  Capus  et  al.  in  [11],  It  has  also  been  applied  to  the  SAR 
and  ISAR  imaging  problems  of  moving  and  maneuvering  targets  by  various  authors  [21, 
22,  23,  24].  The  detection  of  accelerating  targets  in  clutter  is  another  obvious  application 
that  the  FrFT  can  be  expected  to  perform  well,  and  is  a  topic  for  future  reports. 

In  our  problem  of  Doppler  radar  tomography,  the  time  signal  in  a  CPI,  in  many  cases, 
can  be  described  as  the  sum  of  a  number  of  quasi-LFM  signals,  and  thus  the  FrFT  can  be 
effective  in  its  analysis.  In  particular,  if  the  LFM  signals  can  be  extracted  and  spectrally 
compressed  separately  using  a  bilinear  transform  such  as  the  S-Method,  then  the  problem 
with  spurious  components  with  overlapping  signals  can  be  resolved,  and  thus  improved 
image  resolution  can  be  achieved. 


3.1  A  Definition  of  the  FrFT 


Given  a  time  signal  s(t),  its  FrFT  is  defined  as 

/OO 

K(a,u,t)  s(t)  dt,  (24) 

-OO 

where  the  kernel  function  is 

(^^)1/2exp{i^(u2  +  t2)}exp{-i^},  a/nvr, 

K(a,u,t)  =  ^  S(t  —  u),  a  =  n2-K,  (25) 

6(t  +  u),  a  +  Tr  =  n2ir, 


and  where  a  is  called  the  order  of  the  transform,  which  can  also  be  interpreted  as  an  angle 
of  rotation  in  the  time-frequency  plane. 

For  a  =  0  (i.e  for  n  =  0  in  Equation  25),  we  have  Sq(u)  =  s(u )  which  is  the  time- 
domain  signal  itself.  For  a  =  ir/2,  we  have  cot  a  =  0  and  sin  a  =  1,  and  5’T/2('U)  reduces 
to  the  classic  Fourier  transform  S(u)  =  S(f)  of  s(t).  It  is  insightful  to  re-interpret  the 
classic  FT  as  a  rotation  of  the  time  signal  through  an  angle  of  ir/2  in  the  time-frequency 
plane.  A  rotation  through  some  other  angle  gives  a  generalised  or  ‘fractional  frequency’ 
spectrum  of  the  signal  at  that  angle,  with  the  fractional  frequency  denoted  by  u1. 

The  inverse  FrFT  is  simply  the  rotation  of  the  Sa(u)  through  an  angle  of  —a,  such 
that 

/OO 

K(—a,u,t)Sa(u)du.  (26) 

-OO 

where  SQ(u)  =  ^ra{s(t)}. 

xFor  this  reason,  the  same  notation  T  is  used  for  both  the  classic  Fourier  transform  and  the  fractional 
Fourier  transform,  with  a  subscript  a  explicitly  denoting  the  ‘fractional  form’. 
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3.1.1  Properties  of  the  FrFT 

Some  properties  of  the  FrFT  of  order  a,  which  are  most  relevant  to  the  present  work,  are 
summarised  below: 


s(u)  =  JF0{s(t)}  (27) 

The  time-signal  is  its  own  FrFT  at  order  zero. 

2.  Rotation  to  the  classic  Fourier  transform 

S(f)  =  Fn/2  {«(*)}  =  (28) 

The  classic  Fourier  transform  is  the  FrFT  at  order  n/2. 

3.  Cumulative  rotation.  If  a  =  a\  +  «2>  then 

Sa(u)  =  FaAFai{s(t)}}  =  Fai{Fa2{s{t)}}-  (29) 

Successive  rotations  add  up. 

3.1.2  Computation  of  the  FrFT 

We  have  used  a  computational  algorithm  proposed  by  Ozaktas  et  al.  as  described  in  [12] 
for  the  computation  of  Sa(u )  in  (24),  which  can  be  summarised  as  follows. 

1.  First,  the  kernel  function  K(a,u,t )  in  (25),  for  the  cases  of  a  /  rm,  is  re-expressed 
in  a  form  which  is  numerically  more  stable: 

K(a,u,  t)  =  TQexp  j—  hr  u2  tan(  — )  —  f3(u  —  t.)2  +  t2  tan(  — )  |,  (30) 

where  Aa  =  [(1  —  zcot  cc)/27r]1//2,  f3  =  1/sina,  and  u  is  now  in  units  of  Hertz  if 
a  =  7t/2. 

2.  Set  a  value  for  a  and  initialise  (and  discretise)  the  time  vector  t  6  [0 ,tcpi\-  ol  can 
be  restricted  to  an  interval  of  values  corresponding  only  to  practically  realistic  chirp 
rates; 

3.  Multiply  the  original  time  signal  s(t)  with  the  chirp  term  in  the  time  domain, 
exp{— iirt2  tan(a/2)},  ensuring  that  a  sufficient  level  of  interpolation  and  decima¬ 
tion  is  used  to  cover  the  increased  bandwidth  brought  about  by  the  chirp  factor; 

4.  Convolve  the  above  output  with  exp{— m(3t2}] 

5.  Finally,  multiply  the  above  output  with  the  chirp  term  in  the  u  domain  and  the 
complex  factor  Aa. 
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3.2  Application  of  the  FrFT 

3.2.1  Spectral  Composition  of  a  Projection 

Figure  4  depicts  the  ‘chirp  content’  of  a  projection  in  our  Doppler  tomography  problem 
in  which  target  range  is  very  large  compared  to  the  radius  of  the  rotational  motion. 
Theoretically,  each  scattering  center  of  a  rigid  target  traces  out  a  sinusoid  in  the  time- 
frequency  plane  with  the  same  angular  rate  but  generally  different  amplitudes  and  phases 
(which  depend  on  the  relative  position  of  that  center  in  the  target),  and  there  can  be 
many  such  sinusoids.  The  number  of  chirps  in  a  projection  is  equal  to  the  number  of  point 
scatterers  present  in  that  CPI. 

With  the  classic  FT,  a  projection  at  a  particular  point  in  time  can  be  represented  as 
the  FT  spectrum  obtained  from  that  CPI.  The  longer  the  CPI,  the  higher  resolution  can 
be  achieved  in  the  spectral  domain  but  the  frequency  variation  of  the  chirps  also  deviates 
more  from  the  LFM  approximation.  More  importantly,  with  the  classic  FT,  which  is  a 
projection  of  the  signal  onto  the  frequency  axis  (see  Figure  4),  the  chirps  may  significantly 
spread  and  overlap  among  themselves.  With  the  FrFT,  one  or  more  chirps  may  have  their 
spread  much  reduced  in  the  u-domain,  if  a  good  estimate  of  a  is  used  for  the  FrFT,  and 
thus  separation  from  the  rest  of  the  signal  may  be  possible. 


f 


one  projection 

Figure  f:  A  projection  may  contain  many  different  chirps 


To  spectrally  compress  a  chirp  with  chirp  rate  cr  and  its  instantaneous  frequency  in 
the  form  f  =  fo  +  crt ,  the  optimal  transform  order  a  that  should  be  used  is  related  to  cr 
through 

a  =  ^  +  tan_1(T>cr),  (31) 

in  which  the  scale  factor  D  due  to  the  discretisation  of  the  time-frequency  domain  is 
related  to  the  time  and  frequency  resolutions  St  and  Sf  by  [15] 


D  = 


Ns  \ 
PRF2  )  ' 


(32) 
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With  the  angle  of  rotation  expressed  in  (31),  the  signal  would  appear  like  a  ‘pure  tone’ 
along  the  ‘fractional  time’  domain,  or  equivalently  as  an  optimally  compressed  spike  in 
the  ‘fractional  frequency’  domain. 


3.2.2  FrFT  Application 

Since  there  are  many  chirps  of  different  chirp  rates  present  in  a  projection,  a  number  of 
rotations  of  the  time-frequency  plane  may  be  required  for  each  CPI.  Note  that  while  the 
instantaneous  Doppler  frequency  of  a  point  scatterer,  {2u/X)x,  is  determined  by  the  cross 
range  x,  chirp  rate  (derivative  of  frequency)  is  (— 2u>2/\)y,  determined  by  the  down  range 
y  as  seen  in  Figure  5.  The  down  range  extent  of  targets  of  interest  is  finite,  suggesting  a 
finite  range  of  possible  values  over  which  to  search  for  the  ‘best’  or  ‘close  enough’  chirp 
rates. 


Figure  5:  Lines  of  equal  Doppler  frequency  (in  cross-range) 
and  lines  of  equal  chirp  rates  (in  down-range). 


The  relevant  chirp  rate  interval  may  be  sampled  by  a  finite  number  of  values  of  cr.  At 
each  value,  the  fractional  spectrum  is  searched  for  the  most  spectrally  compressed  peak. 

An  illustration  of  the  evolution  of  the  FrFT  for  different  a  values  is  shown  in  Figure  6. 
The  smearing  of  the  spectrum  from  the  two  LFM  signals  and  the  focusing  of  each  spectral 
line  as  the  order  of  the  transform,  a,  is  well  matched  with  the  chirp  rate  of  the  LFM 
signal,  which  occur  at  a  =  84.7°  and  a  =  93.3°,  can  be  seen.  The  traditional  FT  results 
when  a  =  90°,  as  discussed  earlier.  Here,  we  have  used  the  following  system  parameters: 
A  =  0.03  nr,  PRF  =  30  kHz,  rotation  rate  =  300  RPM,  scatterer  1  position  ( x ,  y)  =  (1, 1) 
and  scatterer  2  position  (x,y)  =  (—1,  —1). 


A  search  algorithm  is  required  to  determine  when  the  optimum  focusing  of  a  spectral 
line  occurs,  which  is  then  extracted  by  setting  the  fractional  spectrum  at  all  other  u 
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Fractional  Frequency  (kHz) 


Figure  6:  Evolution  of  the  FrFT  spectrum  of  two  chirp  signals. 


outside  the  focused  peak  to  zero  and  rotated  to  a  =  n/2  (the  physical  frequency  domain). 
The  remaining  portion  of  the  fractional  spectrum  is  rotated  to  the  next  value  of  a  and 
the  process  is  repeated,  until  all  chirp  rate  values  have  been  searched.  A  more  detailed 
description  of  the  algorithm  is  presented  in  Section  5.  All  extracted  chirps  are  then 
spectrally  compressed  using  the  S-Method  described  in  the  next  section. 

The  search  and  extraction  function  is  not  discussed  in  detail  in  this  report  as  it  con¬ 
tinues  to  be  an  area  of  research  and  will  be  addressed  in  future  publications. 


4  The  S-Method  and  its  Application 

In  the  present  algorithm,  the  S-Method  is  used  to  spectrally  compress  a  LFM  signal,  after 
it  has  been  extracted  by  the  FrFT  described  above.  The  S-Method  was  first  proposed  by 
L.  Stankovic  [19]  in  1994  and  has  since  been  applied  to  many  problems  in  micro-Doppler 
signal  analysis  by  T  Thayaparan  et  al.  [6] .  In  the  following,  we  summarise  its  basic  theory 
and  discuss  the  necessary  modification  for  its  application  to  Doppler  radar  tomography. 
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4.1  Description  of  the  S-Method 

The  S-Method  as  proposed  by  Stankovic  is  basically  a  Wigner  distribution2  (WD)  modified 
with  a  reduced  integration  extent. 

While  the  Fourier  transform  is  a  linear  transform  of  a  signal  s(t),  the  Wigner  distri¬ 
bution  is  a  second-order,  or  bilinear  transform,  and  is  defined  (in  the  analog  pseudo  form) 
as 

/OO 

w(t/2)w(-t/2)  s(tk  +  t/2)s*(tk  -  t/ 2)e~l2'Kft  dt.  (33) 

-OO 

It  is  known  that  the  Wigner  distribution  is  ideal  for  a  single  LFM  signal,  in  the  sense  that 
the  spectral  spread  caused  by  the  frequency  modulation  as  seen  in  its  Fourier  transform 
S(f,tk)  is  completely  removed,  and  its  Wigner  distribution  <SVd(/,  Lfc )  is  identical  to  the 
spectrogram  of  a  pure  sinusoid.  For  multiple  LFM  signals,  however,  the  usefulness  of  the 
Wigner  distribution  is  much  reduced  by  the  existence  of  strong  cross  terms. 

The  Wigner  distribution  Swo{f,tk)  can  also  be  expressed  in  terms  of  the  short-time 
Fourier  transform  S(f,tk)  as 

1  r°° 

SwD(f,tk)  =  -  S(f  +  f,tk)S*(f-f',tk)<tf',  (34) 

^  J — OO 

which  can  be  readily  verified  by,  for  example,  directly  substituting  the  definition  of  5(/,  tk) 
as  expressed  in  (13)3.  To  avoid  cross  terms,  at  least  when  the  auto-terms  are  reasonably 
separated  in  frequency,  the  S-Method,  in  its  analog  form,  is  then  defined  as 

1  ff^/2 

SsmU,  tk)  =  ~  S(f  +  /',  tk)  S*(f  -  /',  tk)  df ,  (35) 

77  J-Il/2 

where  fj,  is  a  finite  real  number,  representing  a  finite  frequency  extent  of  integration 
centered  around  the  frequencies  of  the  component  LFM  signals. 

In  discrete  form,  the  S-Method  transform  is  written  as 

L 

SsM(n,k)  =  ^  S(n,k  +  i)  S*(n,k  —  i),  (36) 

i=—L 

where  L  is  an  integer  to  be  chosen.  When  L  =  0,  the  S-Method  reduces  to  the  spectrogram. 
When  L  spans  the  entire  discrete  Fourier  domain,  the  S-Method  is  identical  to  the  Wigner 
distribution.  L  can  be  chosen  to  vary  with  frequency  filter  index  k,  and  since  the  aim 
is  to  optimise  the  focusing  (or  spectral  compression)  of  the  auto-terms  and  to  avoid  the 
cross  terms,  L  needs  to  be  determined  adaptively  for  each  value  of  k.  In  the  present  work, 
the  criterion  used  is  that,  for  a  ‘simple  hump’  in  the  STFT,  the  sum  in  (36)  stops  when 
either  | S(n,k  +  f)|  or  \S(n,  k  —  i)|  is  less  than  a  certain  percentage  of  the  nearest  local 
maximum.  For  multiple  humps,  we  have  used  an  improved  criterion  (not  described  here) 
which  produces  an  L-profile  such  as  in  Figure  8. 

2 which  is  also  known  as  ”  ‘ Wigner-Ville  Distribution’” 

3 Assuming  the  integration  extent  in  (13)  is  from  — oo  to  +oo  and  using  the  property  of  the  Dirac  5 
function. 
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Frequency  Filters 


Figure  1:  Demonstration  of  the  compression  performance  of  the  S- 
Method  technique,  a)  S-Method  and  FFT  comparison,  b)  the  L  profile. 


More  importantly,  to  obtain  a  projection  for  tomographic  imaging  and  preserve  its  first- 
order  dimensionality,  we  also  make  a  simple  modification  of  the  ‘conventional’  S-Method 
by  taking  the  square  root  of  the  above  sum: 


SSM(n,k ) 


L 


1/2 


S(n,  k  +  i)  S*(n,  k  —  i) 

-i=-L 


(37) 


This  modified  S-Method  may  appropriately  be  called  the  tomographic  S-Method. 


4.1.1  Example:  Implementation  of  the  S-Method 

An  example  of  a  smeared  spectrum,  the  ‘spectral  compression  effect’  of  the  S-Method  and 
the  L  profile  is  shown  in  Figure  7.  Note  in  this  figure  that  the  L  profile  is  necessarily 
larger  for  a  more  spread  LFM  signal  (on  the  right)  than  for  the  less  spread  signal  (on  the 
left).  When  a  signal  is  already  strongly  peaked,  no  L  profile  is  necessary. 

4.1.2  Limitation  of  the  S-Method 
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The  S-Method  works  well  for  the  above  example,  when  the  two  chirp  signals  are  well 
separated  in  frequency.  When  the  two  chirp  signals  strongly  overlap,  which  is  certainly 
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Figure  8:  The  compression  performance  of  the  S-Method  technique  for 
multiple,  closely  spaced  terms,  a)  S-Method  and  FFT  comparison  show¬ 
ing  unwanted  cross-terms,  (b)  the  L  profile. 


the  case  when  two  scattering  centers  roughly  line  up  in  the  radar  line-of-sight,  spurious 
frequency  components  arising  from  cross  terms,  which  become  a  nuisance  leading  to  a 
degradation  in  the  final  tomographic  image.  This  undesirable  effect  is  illustrated  in  Figure 
8.  In  general  it  is  desirable  to  apply  an  extraction  step  using  the  FrFT  prior  to  the 
application  of  the  S-Method. 


5  Algorithm  for  Chirp  Extraction  and 

Compression 

This  section  describes  the  algorithm  proposed  and  employed  in  this  work,  for  improving  the 
resolution  of  a  tomographic  image  using  the  fractional  Fourier  transform  and  S-Method. 
Figure  9  shows  a  flow  chart  of  the  algorithm  which  is  described  in  more  detail  below,  with 
Figure  10  providing  descriptive  support. 

Step  1  :  Collect  time  signal  data  for  a  CPI.  Data  from  the  time  signal  is  collected  for  a 
CPI,  with  an  example  shown  in  Figure  10(a). 

Step  2  :  Calculate  the  chirp  rate  and  a  value.  The  first  chirp  rate  is  calculated  based  on  a 
reasonable  assumption  of  the  objects  down  range  extent  and  the  associated  a  is  calculated 
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*  remaining  signal  after  the  chirp  has  been  removed 

Figure  9:  A  flowchart  to  describe  the  processing  steps  to  apply  the  FrFT  and  S-Method 
algorithms  to  data  to  improve  projection  resolution. 
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Figure  10:  FrFT  and  S-Method  algorithm  description  showning  the  process  steps,  a) 
the  time-domain  signal  for  a  CPI  interval,  b)  FrFT  revealing  a  compressed  chirp  and 
uncompressed  chirps,  c)  extracted  compressed  chirp,  d)  FrFT  chirp  rotated  to  the  Fourier 
domain,  e)  S-Method  applied  to  chirp,  f)  S-Method  applied  to  all  chirps  in  CPI,  g)  resulting 
compressed  time-frequency  ( projections )  plot  for  all  CPI. 


using  Equation  31. 

Step  3:  Compute  FrFT  and  chirp  search.  For  the  calculated  a ,  the  fractional  Fourier 
transform  is  applied  to  the  data.  The  resulting  fractional  spectrum  is  searched  for  focused 
chirps,  or  compressed  spectral  lines.  Figure  10(b)  shows  an  example  of  a  fractional  spec¬ 
trum  with  a  compressed  chirp.  If  the  resulting  spectrum  does  not  contain  a  compressed 
spectral  line  then  return  to  step  2  for  the  next  chirp  rate  value.  (The  actual  algorithm  for 
determining  the  existence  of  a  compressed  spectral  line  is  not  discussed  here  and  will  be 
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discussed  in  a  separate  report.) 

Step  4:  Chirp  extraction.  If  a  focused  line  exists,  the  extent  of  this  line  is  extracted  giving 
a  new  spectrum  as  seen  in  Figure  10(c),  which  continues  to  be  processed  by  steps  5  to  7 
below.  The  remnant  signal  goes  back  to  step  2  with  the  next  chirp  rate  value,  until  all 
chirp  rates  have  been  searched. 

Step  5  :  Rotate  to  the  classic  frequency  domain.  At  the  input  to  this  stage,  the  extracted 
chirp  is  still  in  the  fractional  frequency  domain  and  must  be  rotated  to  the  physical 
frequency  domain.  Its  resulting  spectrum  is  an  uncompressed  chirp  as  shown  in  Figure 
10(d). 

Step  6:  Apply  the  S-Method.  The  S-Method  can  now  be  applied  to  the  single  chirp 
spectrum  to  refocus  the  spectral  content.  The  ability  of  the  S-Method  to  properly  focus 
the  chirp  is  improved  by  the  isolation  of  chirps.  The  resulting  S-Method  compressed 
spectrum  is  shown  in  Figure  10(e). 

Step  7:  Compute  the  Projections.  The  isolation  and  extraction  of  chirps  in  the  signal 
are  applied  for  all  chirp  rates  and  the  S-Method  compression  is  applied  to  each.  The  final 
projection  is  a  sum  of  all  compressed  spectra,  as  shown  in  Figure  10(f),  for  a  given  CPI. 

The  processing  steps  described  above  are  repeated  for  all  CPFs  to  give  focussed  pro¬ 
jections  for  all  aspect  angles,  as  exemplified  in  Figure  10(g),  used  later  for  tomographic 
imaging.  Note  that  the  best  estimate  of  the  chirp  rate  and  therefore  of  the  rotation  rate, 
oj,  can  be  determined  by  a  refinement  of  step  2.  However,  this  is  not  the  goal  of  this 
current  work. 


6  Algorithm  Performance  Results 

This  section  aims  to  demonstrate  the  performance  of  the  proposed  image  resolution  im¬ 
provement  algorithm  on  chirped  signals.  The  demonstration  uses  3  types  of  processing. 
The  first  type  uses  an  optimal  number  of  samples,  Ns.  that  gives  no  smearing,  but  rel¬ 
atively  low  resolution,  in  the  traditional  Fourier  transform  based  scheme.  The  second 
type  uses  a  number  of  samples  several  times  greater  than  Ns.  producing  a  smeared  im¬ 
age.  The  third  type  has  the  new  algorithm  applied  to  the  smeared  signal  to  produce  a 
high-resolution  image. 

The  demonstration  is  first  applied  to  rotating  point-scatterers  orientated  in  a  variety 
of  configurations,  and  in  the  presence  of  noise,  and  then  to  real  helicopter  data. 


6.1  Point-scatterer  Model  Target 

6.1.1  Noise-free  Performance 

This  section  will  present  several  examples  to  demonstrate  the  improved  resolution  and 
performance  of  the  algorithm,  in  the  case  of  no  additive  noise  in  the  received  signal.  Table 
1  shows  the  radar  and  object  parameters  used  for  generating  smear-free  and  smeared 
projections. 
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Table  1:  Table  of  radar  and  object  parameters  for  smeared  and 
smear-free  projection  generation. 


Radar  and  Object  Parameters 

Carrier  Frequency  (GHz) 

10 

Pulse  Repetition  Frequency  (kHz) 

35 

Dwell  Time  (s) 

0.3 

Rotation  Rate  (RPM) 

300 

Maximum  Blade  Length  (nr) 

6 

Tip  Doppler  (kHz) 

12.56 

Smear-free  Parameters 

Number  of  samples,  Ns 

57 

TCpi  (nrs) 

1.6 

A 0  (deg) 

3 

Smeared  Parameters 

Number  of  samples,  Ns 

456 

TCPi  (nrs) 

13 

Ad  (deg) 

23 

Smear-free  and  smeared  projections  and  associated  tomographic  images 

A  set  of  4  rotating  point  scatterers  orientated  in  a  square  separated  by  8  nr,  is  shown  in 
Figure  11(a).  The  projections  were  calculated  using  the  value  of  Ns,  given  by  Equation 
19,  leading  to  the  time- frequency  plot  in  Figure  11(b).  Here  we  can  see  the  four  sinusoids 
as  they  rotate  over  time.  Figure  11(c)  shows  the  resulting  tomographic  image,  giving  four 
points  in  the  image.  The  low  value  of  Ns  gives  a  rather  low  resolution  image  but  free  of 
smearing. 

Figure  12(a)  and  Figure  12(b)  show  the  smeared  projections  and  resulting  smeared 
tomographic  image.  The  high  number  of  samples,  8  x  Ns,  used  for  a  projection  gives 
a  finer  image  but  the  longer  integration  time  produces  frequency  smearing  causing  the 
points  in  the  image  to  be  smeared  in  the  tangential  direction. 

Figure  12(c)  and  Figure  12(d)  shows  the  improved  projections  and  resulting  tomo¬ 
graphic  image  after  the  fractional  Fourier  transform  and  S-Method  algorithm  is  applied  to 
the  smeared  projections  to  focus  the  chirps  in  the  signal.  As  can  be  seen,  the  smearing  ef¬ 
fect  is  removed  giving  4  sharp  points  in  the  image.  This  allows  us  to  use  longer  processing 
intervals,  giving  greater  frequency  resolution,  without  smearing  in  the  image. 

To  demonstrate  in  more  detail  the  improvement  using  the  proposed  FrFT  and  S- 
Method  algorithm,  Figure  13  shows  smeared  focused  projections  containing  2  point  scat¬ 
terers  for  three  CPI’s,  especially  when  the  spectra  of  the  point  scatterers  become  more  and 
more  overlapped.  As  can  be  seen,  no  spurious  components  are  created.  Other  illustrative 
examples  are  provided  below. 

Example  1:  Closely-spaced  point  scatterers  near  the  center  of  rotation.  A  set  of  rotating 
point  scatterers  orientated  in  a  square,  starting  at  the  origin  and  separated  by  0.5  nr,  is 
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(c)  Smear- free  tomographic  image. 
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Figure  11:  Projections  and  tomographic  image  using  the  smear-free  Ns. 


Frequency  (kHz)  Frequency  (kHz) 


DSTO-RR-0357 


Cross  Range  (m) 


(a)  Smeared  projections 


(b)  Smeared  tomographic  image 


350  • 

300  - 

250  • 

200  V  3  - 

150  • 

100  - 

50  • 

20  40  60  80  100  120  140 

Time  (secs) 

(c)  FrFT  and  S-Method  improved  pro¬ 
jections 


(d)  FrFT  and  S-Method  improved  to¬ 
mographic  image 


Figure  12:  Smeared  and  FrFT  S-Method  improved  projections  and  tomographic  image. 
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Figure  13:  Smeared  and  FrFT  S-Method  improved  projections  at  three  CPFs,  showing 
the  improvement  in  resolution  and  ability  to  resolve  closely  spaced  scatterers. 


shown  in  Figure  14(a).  The  smear-free  tomographic  image  is  shown  in  Figure  14(b). 
As  expected  we  seen  four  low  resolution  scatterers.  The  smeared  tomographic  images  are 
shown  in  Figure  14(c).  As  can  be  seen,  the  scatterer  at  the  origin,  having  no  radial  extent 
(i.e.  at  the  center  of  rotation),  gives  no  smearing  and  the  overall  quality  of  the  image  is 
poor.  Figure  14(d)  shows  the  FrFT  and  S-Method  improved  tomographic  image,  giving 
four  points  in  the  image.  The  proposed  algorithm  can  easily  focus  4  closely  spaced  point 
scatterers. 

Example  2:  Closely-spaced  point  scatterers  offset  from  the  center  of  rotation  by  4  m.  We 
offset  the  square  by  a  radius  of  4  nr  as  shown  in  Figure  15.  The  low  resolution  in  the  snrear- 
free  image  gives  no  information  about  the  number  of  scatterers  or  the  scatterer  position  as 
seen  in  Figure  15(b).  The  smeared  image  caused  by  the  large  radial  extent,  causes  extensive 
interference  of  the  closely  spaced  scatterers  giving  a  smeared  tomographic  image  that  also 
doesn’t  resemble  the  scatterer  positions.  The  improved  image  after  processing,  distinctly 
shows  the  scatterers  with  high  enough  resolution  to  confidently  resolve  the  scatterers 
positions. 

Example  3:  Linearly  spaced  rotating  point  scatterers  orientated  tangentially  to  the  direc¬ 
tion  of  rotation.  In  this  example,  a  set  of  4  point  scatterers  in  a  linear  pattern  orientated 
tangentially  to  the  path  of  rotation,  separated  by  0.2  m,  as  shown  in  Figure  16(a).  Here 
we  position  the  scatterers  such  that  the  separation  cannot  be  resolved  using  the  snrear-free 
projections,  as  seen  in  Figure  16(b),  and  smeared  projections  give  an  image  with  heavily 


24 


DSTO-RR-0357 


□ 

center  of  rotation  — ' 

05 

04 

03 

02 

01 


(a)  Scatterers  orientation 


(b)  Smear-free  tomographic  image 


(c)  Smeared  tomographic  image  (d)  FrFT  and  S-Method  improved  to¬ 

mographic  image 


Figure  14:  Smeared  and  improved  tomographic  images  for  closely  spaced  rectangularly 
orientated  rotating  point  scatterers 


overlapping  points,  which  is  indeed  the  case  as  seen  in  Figure  16(c).  The  improved  image, 
shown  in  Figure  16(d),  demonstrates  the  excellent  performance  of  the  proposed  algorithm 
even  for  scatterers  position  0.2  m  apart  and  in  a  tangential  configuration  that  would  be 
most  smeared  in  conventional  FT-based  processing. 


Example  4:  Point  scatterers  orientated  in  an  ellipse  offset  from  the  center  of  rotation. 
A  set  of  10  point  scatterers  orientated  in  an  ellipse  offset  from  the  center  of  rotation  is 
shown  in  Figure  17(a).  The  smear-free  image  seen  in  Figure  17(b),  shows  the  position 
of  the  scatterers  but  at  low  resolution.  Using  the  longer  integration  time  (Figure  17(c)) 
to  improve  resolution,  the  smeared  image  shows  how  smearing  increases  as  the  scatterer 
radial  extent  increases,  with  closely  spaced  scatterers  further  away  from  the  center  of 
rotation  overlapping.  The  FrFT  and  S-Method  processed  projections  gives  an  improved 
tomographic  image  (Figure  17(d))  showing  the  scatterers  positions  with  higher  resolution. 
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Figure  15:  Smeared  and  FrFT  S-Method  improved  tomographic  images  for  closely  spaced 
square,  offset  by  4m 


6.1.2  Algorithm  Performance  in  the  Presence  of  Noise 

In  this  section,  the  FrFT  and  S-Method  algorithm  is  tested  in  the  presence  of  noise.  Noise 
samples  are  added  to  the  time-domain  data  samples  and  are  modelled  as 

s(t)  =  s0(t)  +  n(t)  (38) 

where  so(t)  is  the  noise-free  signal  and  n{t)  is  Gaussian  i.i.d.  noise  samples  quantified  using 
the  signal-to-noise  ratio  (SNR). 

Using  the  scatterer  orientation  shown  in  Figure  11(a),  the  processed  tomographic  image 
is  presented  for  varying  SNR.  Figure  18(a)  and  18(b)  show  the  smeared  and  processed 
tomographic  images  respectively  for  a  SNR  =  OdB.  The  amount  of  detail  in  the  smeared 
image  is  hard  to  see  with  the  naked  eye  with  the  smearing  reducing  the  dynamic  range  of 
the  image.  Its  almost  impossible  to  determine  the  position  of  the  four  scatterers.  However, 
the  scatterers  are  easily  seen  in  the  processed  image  for  this  level  of  noise.  This  can  be 
explained  by  the  greater  integration  time  (8  x  Ns)  used  with  the  proposed  algorithm. 
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(a)  Scatterers  orientation 


(b)  Smear-free  tomographic  image 


(c)  Smeared  tomographic  image  (d)  FrFT  and  S-Method  improved  to¬ 

mographic  image 


Figure  16:  Smeared  and  FrFT  S-Method  improved  tomographic  images  for  closely  spaced 
linear  orientated  scatterers  perpendicular  to  center  of  rotation 


Figure  18(c)  shows  the  processed  tomographic  image  for  a  SNR  =  — 4dB.  The  smeared 
image  shows  no  detail  and  the  scatterers  are  difficult  to  see.  The  image  is  starting  to 
degrade  with  scatterer  position  becoming  less  defined.  Figure  18(d)  shows  the  processed 
tomographic  image  for  a  SNR,  =  — 6dB.  Here  we  see  a  degradation  of  image  quality  and 
total  destruction  of  scatterer  information. 

The  FrFT  and  S-Method  algorithm  performs  well  in  de-chirping  the  projections  and 
giving  a  high  resolution  image.  The  algorithm,  however,  relies  on  determining  when  a  chirp 
is  compressed,  which  is  Step  3  in  Section  5  of  the  algorithm  description.  Further  work  in 
this  area  is  continuing  to  maximise  the  robustness  and  performance  of  the  algorithm. 

Also,  this  analysis  used  an  integration  time  such  that  the  number  of  samples  used  is  8 
times  the  smear- free  Ns.  It  can  be  argued  that  by  increasing  this  integration  time,  say  to 
16  x  Ns,  further  performance  in  noise  can  be  achieved.  However,  on  implementing  this, 
no  improvement  in  image  quality  was  seen,  with  signals  with  a  SNR  =  —  6dB  still  giving 
a  degraded  image.  The  reason  for  this  is  the  deviation  from  the  LFM  assumption  as  we 
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(a)  Scatterers  orientation 


(b)  Smear-free  tomographic  image 


(c)  Smeared  tomographic  image 


(d)  FrFT  and  S-Method  improved  to¬ 
mographic  image 


Figure  17:  Smeared  and  FrFT  S-Method  improved  tomographic  images  for  point  scatter¬ 
ers  orientated  in  an  ellipse  offset  from  the  center  of  rotation. 


increase  the  integration  time  and  move  into  the  more  sinusoidal  FM  part  of  the  chirp. 


6.2  Real  Helicopter  Data 

In  this  section,  we  apply  the  proposed  algorithm  to  real  data  of  a  couple  of  helicopters, 
despite  the  fact  that  a  robust  algorithm  for  detailing  and  tracking  focused  ‘humps’  of 
scatterers  is  not  yet  available.  The  radar  return  from  a  helicopter,  apart  from  that  of  the 
helicopter  body  itself,  consists  of  LFM  chirp  signals  from  the  rotating  rotor  blades  with 
the  return  from  a  blade  body  dominating  when  orthogonal  to  the  radar  line-of-sight.  In 
fact,  due  to  the  shape  of  the  rotor  blade  tips,  the  return  from  the  blade  tip  is  consistent 
for  a  wide  range  of  aspect  angles  and  is  indicative  of  a  rotating  point  scatterer  described 
earlier.  Hence  the  techniques  and  FrFT  S-Method  algorithm  described  can  be  applicable 
to  real  helicopter  data  in  an  attempt  to  image  the  blade  tips.  The  two  helicopter  data 
sets  used  for  this  investigation  are  shown  in  Table  2. 
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(a)  Smeared  tomographic  image  with 
SNR  =  OdB 


(c)  FrFT  and  S-Method  improved  to¬ 
mographic  image  with  SNR  =  — 4dB 


(b)  FrFT  and  S-Method  improved  to¬ 
mographic  image  with  SNR  =  OdB 


(d)  FrFT  and  S-Method  improved  to¬ 
mographic  image  with  SNR  =  — 6dB 


Figure  18:  Processed  tomographic  images  of  the  four  scatterer  orientation  for  three  noise 
levels. 


Figure  19  shows  a  3D  view  of  the  tomographic  image  using  the  smear-free  projections 
of  aircraft  A.  The  image  shows  the  peaks  of  the  four  rotor  blades  tips,  but  with  low 
resolution.  However,  Figure  20  is  the  same  view  of  the  same  aircraft  with  FrFT  and 
S-Method  processing  using  smeared  projections.  This  figure  clearly  shows  four  peaks 
corresponding  to  the  four  blade  tip  scatterers.  The  level  of  noise  on  the  signal  is  low 
enough  for  the  chirp  search  and  extraction  algorithm  to  function  properly  and  hence  gives 
a  reasonably  good  tomographic  image. 

Similarly,  Figure  21  shows  a  3D  view  of  the  tomographic  image  using  smear- free  projec¬ 
tions  of  aircraft  B.  Again,  the  low  resolution  and  noise  level  gives  a  less  defined  tomographic 
image.  The  processed  smeared  projections  using  the  FrFT  and  S-Method  algorithm,  for 
the  same  dataset,  exaggerates  the  two  blade  tip  scatterers,  as  shown  in  Figure  22.  Al¬ 
though  the  image  is  less  defined  when  compared  with  the  examples  in  Section  6.1,  the  tip 
scatterers  can  be  easily  seen  and  the  existence  of  two  blade  tips  can  be  verified. 
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Table  2:  Radar  and  helicopter  specific  parameters  of  real  data 
used  to  test  the  FrFT  and  S-Method  algorithm. 


Parameters 

Aircn 

A 

ift  ID 

B 

Pulse  Repetition  Freq.,  FRF 

16.7  kHz 

28.9  kHz 

Carrier  Frequency,  fc 

5.5  GHz 

9.7  GHz 

Rotation  Rate,  u> 

230  RPM 

350  RPM 

Number  of  Blades 

4 

2 

Overlapping  Factor,  F 

0.8 

0.8 

The  FrFT  and  S-Method  algorithm  has  proven  to  perform  well  on  real  helicopter  data 
with  the  tip  scatterers  easily  recognisable  in  the  tomographic  image.  In  less  noisy  data  and 
with  an  improvement  of  chirp  search  and  extraction  algorithms,  the  FrFT  and  S-Method 
technique  could  prove  useful  for  imaging  helicopter  rotor  blades  and  other  rotating  objects. 


7  Conclusion 

We  have  presented  a  new  method  of  improving  the  resolution  of  a  tomographic  image  of 
a  rotating  object  by  effectively  increasing  coherent  integration  time  and  using  recently 
developed  signal  processing  techniques. 

The  imaging  approach  is  Doppler  radar  tomography,  with  the  Doppler-processed  cross¬ 
range  profiles  as  projections  in  the  tomographic  processing.  We  have  briefly  reviewed 
the  basics  of  this  processing  technique  and  determined  its  performance  limit,  in  the  case 
when  only  the  Fourier  transform  is  employed.  The  recently  developed  techniques  are 
the  fractional  Fourier  transform  and  the  S-Method,  here  for  the  first  time,  applied  to 
Doppler  radar  tomography.  The  usefulness  of  these  techniques  lie  in  their  ability  to  identify 
and  spectrally  compress  chirped  signals  which  result  from  scattering  centers  of  a  rotating 
object.  The  projections  are  ‘cleaner’  and  sharper,  resulting  in  tomographic  images  with 
better  resolution. 

Separately,  the  role  of  the  fractional  Fourier  transform  is  to  identify  the  different  chirp 
signals  present  in  a  CPI,  based  on  the  property  that,  at  a  particular  value  of  the  order 
of  the  transform,  only  those  chirps  with  chirp  rates  matching  the  transform  order  can  be 
‘de-chirped’.  An  extraction  algorithm  is  then  used  to  separate  them  from  the  other  chirps 
in  the  signal.  The  role  of  the  S-Method  is  to  spectrally  compress  the  extracted  chirps  in 
the  classic  frequency  domain,  without  the  undesirable  effects  of  cross  terms. 

We  have  demonstrated  the  performance  of  this  new  method  on  a  variety  of  rotating 
point  scatterer  examples,  in  varying  orientations,  in  terms  of  improved  resolution  and  in 
the  presence  of  noise.  We  have  also  demonstrated  the  applicability  of  this  method  on  real 
helicopter  data.  Despite  the  fact  the  current  form  of  the  method  is  not  yet  optimised  for 
the  complexity  of  blade  flashes  from  helicopter  rotor  blades,  the  result  is  very  promising, 
with  the  improved  tomographic  images  showing  focused  blade  tips,  compared  to  smeared 
images  generated  using  Fourier  transform  based  techniques. 
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Figure  19:  3D  smear-free  tomographic  image  of  aircraft  A. 


(m) 


Figure  20:  3D  FrFT  and  S-Method  processed  tomographic  image  of  aircraft  A. 


31 


Magnitude 


DSTO-RR-0357 


Figure  21:  3D  smear-free  tomographic  image  of  aircraft  B. 


Figure  22:  3D  FrFT  and  S-Method  processed  tomographic  image  of  aircraft  B. 
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Other  potential  applications  include  the  imaging  of  helicopter  tail  rotor  blades,  pro¬ 
pellers  on  propeller  driven  aircraft,  and  slower  rotating  objects  such  as  rotating  antenna 
structures  on  ships.  Work  is  currently  underway  to  maximise  the  robustness  of  the  pro¬ 
posed  algorithm  for  many  real  objects  of  interest. 
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